Title of the Invention: 

MYOELECTRIC FEATURE -PAT TERN CLASSIFICATION 
METHOD AND APPARATUS 

BACKGROUND OF THE INVENTION 
Field of the Invention: 
[0001] 

The present invention relates to a myoelectric 
feature-pattern classification method and apparatus in a 
method of interfacing muscle action potential (myoelectric 
pattern) . 

Description of the Prior Art: 
[0002] 

Figure 2 is a drawing used to illustrate a prior-art 
apparatus that manipulates a control target by extracting a 
feature pattern from a myoelectric pattern and classifying 
the extracted feature pattern. In the drawing, CD denotes 
myoelectric patterns, © a surface electrode group, © 
amplification and smoothing apparatuses, © a feature-pattern 
extraction apparatus, © a pattern classifier, and © control 
target such as a motor, a robot, a device for the disabled, 
a rehabilitation device, a myoelectric arm prosthesis, a game, 
and so forth. 
[0003] 

As shown in the drawing, a myoelectric pattern © that 
is an action potential generated by the coordinated action 
of a plurality of muscles, is measured by a set or a plurality 
of sets of surface electrode groups © on a skin surface. What 
is measured is the sum of the action potentials generated by 
the plurality of muscles. 
[0004] 

Next, the measured potential is subjected to 
amplification and smoothing by the amplification and smoothing 
apparatuses ©. The feature-pattern extraction apparatus © 
extracts a feature pattern from the amplified, smoothed signal . 
The pattern classifier © classifies feature patterns from the 



extracted signal, and generates control signals to control 
the control target ©. A real value filter, such as a neural 
network or the like, or a logical value filter, such as a logic 
circuit or the like, can be used for the pattern classifier. 
[0005] 

This type of prior-art apparatus requires, as the 
feature-pattern extraction apparatus , a high-level arithmetic 
processing device {a high- specif i cat ion CPU or DSP, or a special 
LSI) for FFT calculations and the like and for solving inverse 
mapping problems, which has been a problem standing in the 
way of reducing the size and cost of the apparatus. This has 
prevented the apparatus coming into widespread use as a 
myoelectric pattern interface. 
[0006] 

Sometimes a bias in myoelectric pattern distributions 
can make it difficult to achieve classification with good 
accuracy. It is possible that the bias in the distribution 
is caused by the distance relationship between the measurement 
electrodes and contracted muscle. That is, because a 
myoelectric pattern at the contraction of a muscle that is 
far from the electrodes is attenuated as it passes through 
the living body tissue, it distributes in low-value regions. 
In contrast to this, in the case of the contraction of a muscle 
that is close to the electrodes, there is little attenuation 
from the propagation in the body, so it distributes in high-value 
regions . 
[0007] 

The present invention is proposed to resolve the above 
problems, and has as an object to provide a myoelectric 
feature-pattern classification method and apparatus that is 
smaller and cheaper, promoting the wider use of a myoelectric 
interface apparatus. 
[0008] 

Another object of the present invention is to provide 
a myoelectric feature-pattern classification method and 
apparatus that, in cases where the myoelectric pattern 



distribution is biased, can reduce the distribution bias and 
increase the pattern classification accuracy. 

SUMMARY OF THE INVENTION 

[0009] 

The myoelectric pattern classification method of the 
present invention comprises using logarithmic transformation 
processing to extract a feature pattern from a myoelectric 
pattern that is a muscle action potential, classifying the 
extracted feature pattern and generating an output control 
signal . 
[0010] 

Also, the myoelectric pattern classification apparatus 
of the present invention comprises a feature-pattern 
extraction apparatus that uses a logarithmic transformation 
apparatus to extract a feature pattern from a myoelectric 
pattern that is a muscle action potential, and a pattern 
classifier that classifies the extracted feature pattern and 
generates an output control signal. 
[0011] 

The distance relationship between a contracted muscle 
and the electrodes can give rise to distributions that are 
biased towards low-value regions and high-value regions. In 
this case, it is possible to decrease the distribution bias 
and raise pattern classification accuracy by using logarithmic 
transformation to transform the patterns to increase the 
low-value region resolution and decrease the high-value region 
resolution. 
[0012] 

The above logarithmic transformation can be realized by 
means of an analogue filter, or by software on a CPU or 
microprocessor, or by a lookup-table or the like. This makes 
it possible to reduce the size and cost, thereby promoting 
the wider use of a myoelectric interface apparatus. 



[0013] 

Other objects and features of the invention will be more 
apparent from the following detailed description of the 
invention based on the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] 

Figure 1 is an explanatory drawing of the configuration 
of a myoelectric feature-pattern classification apparatus 
according to the present invention. 
[0015] 

Figure 2 is an explanatory drawing of the configuration 
of a myoelectric feature-pattern classification apparatus 
according to the prior art. 
[0016] 

Figure 3 is a general block diagram of an apparatus for 
measuring myoelectricity . 
[0017] 

Figure 4(a) is a graph showing the state of a muscle 
contracted for three seconds after being relaxed for two 
seconds . 
[0018] 

Figure 4 (b) is a graph showing an example of a myoelectric 
signal during, in the state of Figure 4 (a) , wrist flexion. 
[0019] 

Figure 4 (c) is a graph showing an example of a myoelectric 
signal that is the myoelectric signal of Figure 4 (b) that has 
been subjected to smoothing processing. 
[0020] 

Figure 5 is a drawing for explaining the feature-pattern 
extraction method. 
[0021] 

Figure 6 is a diagramof the distribution of the myoelectric 
feature-pattern extraction method patterns. 



[0022] 

Figure 7 is a diagramof the distribution of the myoelectric 
feature patterns that are the feature patterns of Figure 6 
that have been subjected to logarithmic transformation. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0023] 

Figure 1 is an explanatory drawing of the configuration 
of an apparatus of the present invention that extracts a feature 
pattern from myoelectric patterns and controls the motors used 
to operate an arm prosthesis as the control target ©. In the 
drawing, © denotes myoelectric patterns, © a surf ace electrode 
group, (D amplification and smoothing apparatuses, © a 
feature-pattern extraction apparatus that uses a logarithmic 
transformation apparatus, © a pattern classifier, and © the 
control target. 
[0024] 

As shown in the drawing, a myoelectric pattern © that 
is an action potential generated by the coordinated action 
of a plurality of muscles is measured by one or a plurality 
of surface electrode groups @ on a skin surface. What is 
measured is the sum of the action potentials generated by the 
plurality of muscles. Next, the measured potential is 
subjected to amplification and smoothing by the amplification 
and smoothing apparatuses ©. 
[0025] 

The present invention uses a logarithmic transformation 
apparatus as the feature-pattern extraction apparatus ® . The 
logarithmic transformation apparatus extracts feature 
patterns by performing logarithmic transformation of the 
myoelectric patterns . This logarithmic transformation can be 
realized by means of an analogue filter, or by software on 
a CPU or microprocessor, or by a lookup-table or the like. 
The pattern classifier (D classifies feature patterns from the 
extracted feature patterns, and generates signals to control 
the control target ©. A real value filter, such as a neural 



network or the like, or a logical value filter, such, as a logic 
circuit or the like, can be used for the pattern classifier. 
[0026] 

In order to show the effect of the present invention, 
there will now be explained an example in which classification 
accuracy is improvedby using the method of the present invention 
in the action-decisions of a myoelectric arm prosthesis. 
[0027] 

First, myoelectricity measurement will be described, 
using Figure 3. Because the myoelectric signal is weak, in 
the order of several |J.V to several mV, in the measurement thereof, 
it is generally necessary to subject a signal measured using 
two electrodes to differential amplification to reduce the 
effect of noise generated by power supplies or other electric 
devices. In the illustrated example, the myoelectricity 
measurement is made using eight (four sets) electrodes 12 
applied to four locations around the forearm 11 of a myoelectric 
arm prosthesis. Reference numeral 13 denotes a reference 
electrode . 
[0028] 

The weak myoelectric signals, measured by the four sets 
of electrodes 12 are sent, via a sensor box 14, to an amplifier 
15 where they are differentially amplified and are then sent 
to anA/D converter 16. At the A/D converter 16, the amplified 
myoelectric signals are separated by 1000 Hz and converted 
into 12-bit digital signals and stored on a notebook type 
personal computer (PC) 17 connected by a SCSI interface. In 
the drawing, reference numeral 18 denotes a power supply and 
numeral 19 a ground. In this measurement, also, an analogue 
filter was used that blocks the commercial power supply 
frequency (50 Hz) in order to eliminate noise generated by 
the commercial power supply. 
[0029] 

In the following explanation, equation 1 expresses a set 
of signals measured at time t. 



[0030] 

[Equation 1] 



f(t) - (fi(t)M)Js(t)M)) 

[0031] 

Myoelectricity measurements were carried out while the 
muscles were contracted as if to perform six actions (forearm 
pronation, f orearmsupination, wrist flexion, wrist extension, 
hand closing and hand opening) . That is, myoelectricity 
measurements were carried out while performing the muscle 
contractions shown in Figure 4(a) (contracting the muscles 
after a two-second relaxation period and maintaining the 
contraction for three seconds) twenty times for each of the 
six actions (for a total of 120 times) . In the case of forearm 
amputees, the myoelectricity measurements were carried out 
when the muscles were contracted based on a pre-amputation 
action image. 
[0032] 

Figure 4 (b) shows an example of myoelectricity measured 
on the surface of the skin over a flexed forearm muscle group 
during wrist flexion. The measured myoelectric signals were 
subjected to smoothing processing to remove high-frequency 
components. Such smoothing processing is usually done with 
an analogue low-pass filter. However, in the case of these 
measurements, storage on the PC was effected without using 
a low-pass filter to avoid loss of information, so the smoothing 
processing was carried out on the PC. Specifically, a smoothed 
rectified value at time t is obtained by the calculation shown 
in equation 2. Figure 4(c) shows an example of a smoothed 
rectified myoelectric signal. 
[0033] 

[Equation 2] 



UK) 1(H) 100 100 

ait) = (£ I/, (t - «•)!.£ \Mt - 01. £ \hit - 01, £ \Mt - 01) 



[0034] 

Myoelectric feature patterns used in action-decisions 
of an arm prosthesis are usually extracted using (1) a method 
in which the extraction is made from the myoelectricitymeasured 
at the time of the initiation of the muscle contraction, or 

(2) a method in which the extraction is made from the 
myoelectricity measured during a state in which muscle 
contraction is maintained (hereinafter called the steady 
state) . In the case of method (1) , there is the advantage that 
the delay from the initiation of muscle contraction until the 
armprosthesis starts anactionis short. However, withrespect 
to the action classification rate, method (2) is known to be 
better than method (1), so here, the feature patterns are 
extracted using steady-state myoelectricity . 

[0035] 

The feature patterns are extracted using smoothed, 
rectified sready-state myoelectric values . This smoothing is 
also usually carried out using an analogue low-pass filter, 
but here, the average value of the values obtained on the PC 
by equation 2 during a period of one second is calculated and 
used as the feature pattern. 
[0036] 

In this feature-pattern extraction, in order to make it 
less tiring for the subject during the myoelectricity 
measurements, ten averaged values are obtained from one muscle 
contraction. That is, as shown in Figure 5, averaged values 
per second, each offset 100 milliseconds starting from one 
second after action initiation, are obtained using the 
calculation shown in equation 3. Here, i is the time at which 
an action is initiated, that is, the time at which the value 
obtained by equation 4 exceeds a pre-set threshold value. 
[0037] 

[Equation 3] 



[0038] 

[Equation 4] 



G(t) = g x (t) + ff2(t) + g :i (t) + g 4 (t) 
[0039] 

However, inpractice, depending on how a muscle contracts, 
there are cases in which the value of equation 4 does not exceed 
the threshold value. Therefore, not all the twenty 
contractions per action can be used for the calculation of 
equation 3. Therefore, as training patterns for training the 
pattern classifier, equation 3 values are used from five of 
the muscle contractions in which the equation 4 value exceeded 
the threshold value (5 [muscle contractions] x 10 [patterns] 
x 6 [actions] = 300 [patterns]). The equation 3 values 
extracted from another five muscle contractions are used as 
test patterns for evaluating the pattern classifier. In the 
following explanation, prepared training patterns are 
expressedby equation5 andprepared test patterns are expressed 
by equation 6. 
[0040] 

[Equation 5] 

Tr org (n) = (Tr orsJ (n), Tr org , 2 (n), Tr orgJ (n), Tr orgA (n)) 
(n = 0,1,-, 299) 

[0041] 

[Equation 6] 

Te org (n) = (Te ors>1 (n), Te orgi2 (n), Te orgJ (n), Te orgJ (n)) 
(n = 0,1, -,299) 

[0042] 

Here, an investigation relating to the characteristics 
of myoelectric-pattern distributions will be followed by a 
discussion of a method that utilizes those characteristics 
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to improve classification accuracy. To investigate the 
characteristics of myoelectric-pattern distributions, Figure 
6 shows a plot, for each action, in which the first components 
of the above-described training patterns are X axis values 
and the second components are Y axis values. In this example, 
it can be seen that the characteristics of the 
myoelectric-pattern distributions are broadly divided into 
two . 
[0043] 

The first characteristic is that the myoelectric patterns 
do not distribute uniformly, but distribute with biases . That 
is, they distribute divided into low- value regions and 
high-value regions, with patterns having a distribution biased 
towards low-value regions (supination +, pronation x, flexion 
*) having a small distance between patterns of each action, 
and feature-patterns distributed in high-value regions 
(extension □, hand closing ■, hand opening O) showing a large 
variance. 
[0044] 

This distribution bias may be caused by the distance 
relationship between the measurement electrodes and the 
contracted muscles. That is, because a myoelectric pattern 
at the time of the contraction of a muscle that is far from 
the electrodes is attenuated as it passes through the living 
body tissue, it distributes in low- value regions . In contrast 
to this, in the case of the contraction of a muscle that is 
close to the electrodes, there is little attenuation from the 
propagation in the body, so it distributes in high-value 
regions . 
[0045] 

The following description relates to a method of improving 
classification accuracy by utilizing this characteristic to 
reduce bias in classification-pattern distributions. To do 
this, first, in order to evaluate classification-pattern 
distributions, there is obtained a "within-class variance 
between-class variance rate," which is usually used as an 
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evaluation index. That is, in the usual pattern classification, 
high-accuracy classification is possible when patterns 
belonging to the same group (hereinafter called class) 
distribute closer together and patterns belonging to different 
classes distribute apart. Accordingly, the within-class 
distribution spread (within-class variance: Equation 7) 
between-class distribution spread (between-class variance: 
Equation 8 ) rate can be used as a distribution evaluation index . 
[0046] 

[Equation 7] 

1 6 

^ = 3oo E E - mi)t ( x - m o 

i=i xe Xi 

X' set of patterns belonging to the ith class, 
rrti'. average vector for Xt 
[0047] 

[Equation 8] 

1 6 

= 300 E 50 * (mi - mfirm - m ) 

Xi- set of patterns belonging to the ith class, 
mc average vector for %> 
m: average vector of all patterns 
[0048] 

Specifically, a higher within-class variance 
between-class variance rate, defined by equation 9, enables 
a feature pattern to be evaluated as better. 
[0049] 

[Equation 9] 
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[0050] 

In accordance with these equations, when the within-class 
variance between-class variance rate of the training pattern 
was obtained using equation 5, the average for 13 subjects 
was 0.48. Below, the logarithmic transformations shown next 
were used to carry out classification-pattern transformations 
with the aim of achieving high-accuracy classification of 
myoelectric patterns by improving the within-class variance 
between-class variance rate. 
[0051] 

Here, in order to improve the characteristic of distances 
between patterns for each action being small in low-value 
regions and patterns distributed in high-value regions having 
a large variance, a method that performs logarithmic 
transformation is proposed. With this method, owing to the 
nature of logarithms, in low-value regions resolution is 
increased, and in high-value regions resolution is decreased. 
[0052] 

Accordingly, if, for example, the distribution shown in 
Figure 6 is subjected to logarithmic transformation, as shown 
in Figure 7 , the variance of patterns distributed in high-value 
regions becomes lower, and the distance between patterns 
distributed in low- value regions becomes higher. 
[0053] 

Thus, the calculations shown in equation 10 and equation 
11 were used to perform logarithmic transformation on feature 
patterns. When the within-class variance between-class 
variance rate was obtained using equation 5, the average for 
the 13 subjects was 0.89. That is, relative to the case in 
which logarithmic transformation is not used, it was possible 
to achieve an approximately 1.85-fold improvement, so it can 
be expected that higher-accuracy classification will be 
possible . 



[0054] 

[Equation 10] 

Trio g (n) = (Triogjfa), Tr log , 2 (n), Tr log3 (n), Tr logA (n)) 

= (-log(Tr orgJ (n)), -log(Tr org2 (n)), -Iog(Tr org , 3 (n)), -log(Tr orgA (n))) 
(n = 0,1, -,299) 

[0055] 

[Equation 11] 

Te log (n) = (Te^jfn), Te log2 (n), T eiog , 3 (n), Te logA (n)) 

= (-log(Te orgJ (n)), -log(Te org , 2 (n)), -log(Te org , 3 (n)), -log(Te orgA (n))) 
(n = 0,1, -,299) 

[0056] 

In the following, training patterns obtained with 
equation 5 and equation 10 are each used to train the pattern 
classifier, and test patterns obtained with equation 6 and 
equation 11 are used to evaluate each classifier. In this 
evaluation, the pattern classifier is neural-network-based, 
and the freely distributed NevProp (Nevada back Propagation) 
is the neural network software package used. NevProp is a 
perceptron simulator that can be used extremely easily. 
[0057] 

A feature of NevProp is its use of a method called 
cross-validation to prevent the phenomenon where the pattern 
classifier loses its generalization ability with respect to 
test patterns as a result of excess adaptation to the training 
patterns (hereinafter called overfitting) . In this method, 
the training patterns are divided into two sets, which are 
used in training to enable correct classification of both 
without loss of generalization ability with respect to test 
patterns. NevProp also repeats the cross-validation five 
times to reduce the bias effect during the division of the 
training patterns. 



[0058] 

The structure of the neural network used for the 
classification has four input-layer nodes, eight 
intermediate-layer nodes, and six output-layer nodes. The 
output-layer nodes correspond to the six prosthetic-arm 
actions, and when a threshold value (0.5) is exceeded, the 
corresponding action is selected. 
[0059] 

The result of classification using this neural network 
showed that with the logarithmic transformation, the 
myoelectric-pattern classification ratio improved by an 
average of 2.7% with respect to thirteen subjects, and by a 
maximum of 23.5% (nine subjects). 
["00 60] 

The distance relationship between a contracted muscle 
and electrodes can give rise to myoelectric patterns 
distributed with biases towards low-value regions and 
high-value regions, inwhichcase there is aproblemof decreased 
classification accuracy . Using logarithmic transformation to 
transf ormpatterns makes it possible to reduce the distribution 
bias by increasing the low-value region resolution and 
decreasing the high-value region resolution. This makes it 
possible to increase the pattern classification accuracy and, 
by also expanding the range of possible application, makes 
it possible to promote the wider use of a myoelectric interface 
apparatus. Also, this method can be realized by means of an 
analogue filter, or by software on a CPU or microprocessor, 
or by a lookup-table or the like, making it possible to reduce 
the size and cost. 



